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The effective action describing the gapless Nambu-Goldstone, or Anderson-Bogohubov, mode of 
. a zero-temperature dilute Fermi gas at unitarity is derived up to next-to-leading order in derivatives 

from the microscopic theory. Apart from a next-to-leading order term that is suppressed in the 
BCS limit, the effective action obtained in the strong-coupling unitary limit is proportional to that 
^ , obtained in the weak-coupling BCS limit. 
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I. INTRODUCTION 



The ability to control the interaction between atoms is unique to experiments on optically trapped ultracold atomic 
gases, unmatched by other condensed matter experiments and unthought of even a few years ago. After successfully 
creating a degenerate Fermi sea in a trapped Fermi gas m , Jin and collaborators were able to tune a dilute ultracold 
. Fermi gas to be near a magnetic Feshbach resonance where small changes in magnetic field strength have pronounced 
^ ■ effects on the two-particle scattering length characterizing the interatomic pair potential in vacuum Using this 
O^, handle, the group produced in 2003 for the first time a condensate of pairs of ferniionic atoms outside the weak- 
coupling Bardeen-Cooper-SchriefFer (BCS) limit of loosely bound Cooper pairs Q, and subsequently managed to 
realize the intriguing crossover from the BCS limit to the limit of tightly bound pairs that form a Bose-Einstein 
Q ' condensate [1]. This crossover, which was studied earlier by a number of theorists is accompanied by a smooth 

change in chemical potential from positive in the BCS limit to negative in the Bose-Einstein condensate (BEC) limit. 
(— ( I The unitary, or strong-coupling, limit of the BCS-BEC crossover, which is in the region with positive chemical 



Q ' potential, is of particular interest [12| . This limit marks the threshold of a bound state in vacuum, where the scattering 
O I length diverges and changes from negative on the BCS side to positive on the BEC side. Because the scattering length 
diverges in this limit, the dilute system has no obvious scale parameter, other than the particle number density n. 
J This implies, for example, the remarkable result that the ground-state energy of this strongly interacting system is 
^ ■ proportional to that of a free Fermi gas, or equivalently, a neutral BCS superconductor at the same density. The 
lO ' absence of a coupling constant, on the one hand, precludes a standard perturbative approach, but on the other makes 
feasible an effective field theory approach. This is because, as for critical phenomena characterized by a diverging 
length scale, the absence of an intrinsic scale gives rise to universal behavior. The effective field theory program 
" [ : differs from perturbation theory in that it is carried out not by expanding in a small interaction-related parameter, 
' but by expanding in powers of energy and momentum instead. Using general coordinate and conformal invariance 
, as guiding principle. Son and Wingate [Tsj . in an original paper, proposed the most general effective Lagrangian to 
' leading and next-to-leading orders in a gradient, or momentum, expansion. The conformal transformations involve 
T\ ] arbitrary reparametrizations of time t — t'{t), while the general coordinate transformations are restricted to curved 
^ ■ space (as opposed to spacetime). By keeping the metric gij — 1,2,3) arbitrary in intermediate steps, they were 
able to derive nontrivial results that survive the limit of flat space. At leading order, these new symmetry arguments 
are not more powerful than those based on just Galilei invariance, but they are claimed to be more powerful at the 
' next-to- leading order [isj . 

The importance of Galilei invariance in describing (clean) BCS superconductors at the absolute zero of temperature 
was already stressed in the 1960ies by Kemoklidze and Pitaevskii [IJ], following a suggestion by Nozieres. It has been 
used as a guiding principle for obtaining effective theories of zero-temperature neutral superconductors as well as 
superfluids by others since (T5l - [l8| . Such effective theories are expressed in terms of the Nambu-Goldstone mode 
emerging from the spontaneous breakdown of the global U(l) phase symmetry in such systems. Since this mode is 
gapless, it constitutes the most important degree of freedom at low frequencies and long wave lengths. Its existence 
in neutral BCS superconductors was first pointed out by Anderson fl^ and Bogoliubov |20|. Invariance under global 
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phase transformations implies that a Nambu-Goldstone field is invariably accompanied by at least one derivative. 
Galilei invariance then restricts how gradient terms can appear in combination with terms involving time derivatives. 

In this paper, the effective field theory (EFT) describing a dilute Fermi gas at unitarity and at the absolute zero of 
temperature up to next-to-leading order is derived from the microscopic theory using a derivative expansion method. 
This EFT program was first carried out in Ref. in the weak-coupling BCS limit and in Ref. 22] for a weakly 
interacting Bose gas, and is extended here to the strong-coupling unitary limit, for details see the textbook [23j. Due 
to a vanishing vertex and conspiring contributions, coefficients in the EFT at unitarity can, against common beliefs, be 
computed analytically beyond the Gaussian approximation — albeit approximately. By construction, the approach by 
Son and Wingate [l3[ is limited to an expansion in momentum (or inverse wave length). Since the derivative expansion 
method [H, |2^ we adopt treats time and spatial derivatives on equal footing, it does not face this limitation and 
yields an expansion in both momentum and energy (or frequency). This leads to two additional terms in the EFT at 
next-to- leading order omitted in Ref. [l^. As far as static response functions are concerned, these additional terms 
are immaterial. But they are relevant for the spectrum of the Nambu-Goldstone mode, or phonons, at next-to-leading 
order and modify the result obtained in Ref. [13] . 

The paper is organized as follows. The next section sets the stage and briefly introduces the derivative expansion 
method we adopt. Sections IIIII and IVl treat the effective theory at leading and next-to-leading order, respectively. 
Section HVl discusses the one-dimensional case for which exact results are known. Finally, integrals and vertices needed 
in this study are collected in the Appendix. 

For notational convenience, we adopt a relativistic notation. A spacetime point will be indicated by the four-vector 
X ^ ~ {t, x), /i = 0, 1, 2, 3, while the frequency w and momentum k of a particle will be denoted by fc'^ = (w, k). 
The time derivative dt = d/dt and the gradient V are combined in a single vector 9^ = (9t, V). Indices are raised 
and lowered with the help of the diagonal metric 77^^ = diag(l, —1, —1, —1), so that, for example, = {dt, —V). We 
also write k ■ x = fc^x^ — fc^a;^ for — k ■ x, and use Einstein's summation convention. Natural units h — c = 1 are 
adopted throughout. 



II. DERIVATIVE EXPANSION 

A dilute Fermi gas at the absolute zero of temperature can be modeled by the Lagrangian density 



a=f,l 



/ 9 \ 



where '0t(4.) is an anti-commuting field that describes the fermionic atoms of mass m and spin up (down) , and /i is 
the chemical potential. The true interatomic pair potential, which has typically a repulsive hard core of radius less 
than one nanometer and a weak long-range attractive tail, is approximated by a local, i.e., delta- function potential, 
characterized by a single interaction parameter A. This parameter is related to the (s-wave) scattering length a of the 
true potential, parameterizing two-particle scattering at low energy in vacuum, through 

1 _ V(D/2 - 1) m 

A 47r^/2 ^ > 

in dimensional regularization. For later convenience, we have recorded the result for arbitrary number of space 
dimensions _D, and F denotes the gamma function. For D = 3, Eq. ^ reduces to 1/A = m/47ra. The Lagrangian ([T]) 
is invariant under Galilei transformations, where it is recalled that under a Galilei boost with a constant velocity u, 
the coordinates transform as 

i ^ t' = t, X ^ x' = X - ui (3) 

so that 



d d dt d dx 
dt^ df^ dt'dt ^ dt' 

and ipi^) picks up an extra phase factor 



V = (9t + u • V, V ^ V = V, (4) 



V'(x) ^ ip'ix') = e""(-"-''+2" *) ^(x). (5) 
After a Hubbard-Stratonovich transformation, the zero-temperature partition function, 

Z = I D^^D^-exp (i / d^xc) , (6) 
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can be written in the standard form quadratic in the fermion fields at the expense of additional integrals over auxiliary 
fields A and A^: 



Z = J D*tD* J DA*DA exp ^i^ J d'^x\A{x) 

i/d^.v,tfi9* + VV2m + ^ -A(x) 
J \ -A*{x) idt-V /2m-^ 



X exp 



(7) 



where ^ stands for the two-component Nambu spinor 

*^(j!t), ^t^iv^V'^). (8) 

The fermionic degrees of freedom can now be integrated out exactly with the result 

Z ^ j T)A*T)Ae'^^^'^^\ (9) 
where 5'[A*, A] denotes the one-fermion-loop effective action 



(10) 



with ^(k) — k^/27Ti — /i the kinetic energy of noninteracting fermions measured relative to the chemical potential. 
The trace appearing here is evaluated by using plane waves as a basis 

Trln{/^[p,A(x)]} =tr j d^x J e"=-^ In {T^b, A(a;]} e'"= ■^ (11) 

with "tr" denoting the trace over the discrete indices. We adopt the convention that the momentum operator 
pf^ = {idt, — iV) acts on all the fields to its right, whereas the ordinary derivative 9^ — {dt, —V) acts only on the next 
field to its right. The integral J d'^k stands for the integral over loop momenta fc^ = (^^ik). 

For a static uniform system, where the order parameter A(a;) = A is independent of spacetime, the trace in Eq. (jlip 
reduces to 



Trln[if(p,A)] =tr J d^x J ^ ln[i^(fc,A)] 



(12) 



with the integral over spacetime giving just a volume factor. The frequency integral in Eq. (jlO[) can be evaluated in 
closed form to give for the one-fermion-loop effective action 5[A, A*] = — J dtd'^xVcS with 



Vcff = -i|Ap - / im) - e(k)] , (13) 



the effective energy density and 



E(k) ^ y/eik) + |A|2 (14) 
the single-fermion excitation spectrum of the noninteracting system. The celebrated BCS gap equation, 

1 1 r d3. 1 ^^^^ 



A 2 J (27r)3£;(k)' 

follows by minimizing the effective potential with respect to A. The integral appearing here can be evaluated analyt- 
ically to give [HI 

1 _ 4 hjxo) - XoIe{xo) 
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where Xo = A*/^: ^-nd /s and /g are two integrals which are related to the response functions d'^V^/dfi^ , or d'^Vcs/dA'^ 
and d^Ves/dAdfi, respectively, and which can be expressed in terms of Legendre functions [2y|, see the Appendix. 
The momentum fcp, defined through 

" = ^'^F, (17) 

is used to remove the dimension of the scattering length introduced in Eq. ([2]). This parameter relates the particle 
number density n of the interacting Fermi gas to the Fermi momentum of a free Fermi gas at the same density. The 
density n of the interacting system is obtained from the effective potential (jl3[) by differentiating it with (minus) the 
chemical potential. The result can be put in the form 

- = ^ 273' (18) 

ep [xohiXo) + l6{xo)r 

with ep = fc|,/2m = (37r^n)^/'^/2m the Fermi energy of a free gas. Equations and (IT5|) are valid in the entire 
BCS-BEC crossover, with the BCS limit > 0, a < 0) corresponding to Xo » 1, where 

h{xo) ~ xl/^, h{xo) ~ ln(8xo)/2xy2, (19) 

and the EEC limit (/it < 0, a > 0) corresponding to Xq << —1, where 

h{xo) « 7r/16|xo|3/2, Ie{xo) « 7r/4|xo|i/2. (20) 

At unitarity, where a bound state appears and l/fcpa tends to zero, — XqIq by Eq. (|16p . so that Xo = 0.8604. . ., 
h = 0.8693 . . . , /e = 1.010 . . ., and 

— i —= 0.6864.... (21) 

[{i + xDhf 

For the ratio ^ of the chemical potential and the Fermi energy ep = kp/2m of the free gas, this gives the value 

i^ — =^o— = — ^= 0.5906.... (22) 



"ep [(l+a;2)/6l2/3 



These mean-field values, which were first obtained numerically in Ref. [lO|, should be compared with, for example, the 
estimates A/ep = 0.84(4) and ^ = 0.42(1) obtained through quantum Monte Carlo simulations of systems of about 
60 particles [23]. 

One of the observations of Ref. |l3| is that all the leading order terms in the effective theory describing the Anderson- 
Bogoliubov mode at unitarity can be determined from the expression (jl3p for the static uniform system. This will be 
demonstrated in the following section by explicit calculation. 

To determine the next-to-leading terms, the auxiliary fields A and A^ can no longer be assumed to be constant in 
the formal expression pop . It can then at best be evaluated in a derivative expansion. The method [2^ l25j we adopt 
proceeds as follows. First, the logarithm is expanded in a Taylor series. Each term in the series contains powers of 
the derivative operating on every field appearing to the right. Second, all these operators are shifted to the left 
by repeated use of the identity 

Mx)p''Mx) = (P^ - id^)Mx)Mx), (23) 

where 0i(x) and 4>2ix) are arbitrary fields, and the derivative = {dt, —V) acts by convention only on the next field 
to the right. Next, each term is integrated by parts so that all the p^'s act to the left where only a factor exp(i/c • x) 
is present and yield a factor of k^. In this way, each occurrence of the operator p^ is replaced with an integration 
variable k^. Finally, the exponential function exp(— i/c • x) is moved to the left where it is multiplied with exp(i/c • x) 
to give unity. The momentum integration can now in principle be carried out to yield an effective action written as 
a spacetime integral over a Lagrangian density, S = J d'^x C 

III. LEADING ORDER 

To derive the effective action, we write the complex field A(a::) in terms of a spacetime-dependent amplitude and 
phase as 

A(a;) = [A + (7(2;)]e2'^(^', (24) 



5 



with A denoting a real solution of the gap equation ([T5|) . The functional measure in Eq. ^ must be changed 
accordingly by expressing it in terms of the new fields. Physically, the phase f{x) of the order parameter describes 
the Nambu-Goldstone mode accompanying the spontaneous breakdown of global U(l) phase symmetry. The effective 
action governing these phonons obtains after integrating out the a field. 

The phase can be removed from the order parameter by returning to the expression ([7]) for the partition function 
and decompose the fermion fields as 

^,ix)^e"^'^-'^xAx). (25) 
Instead of the one- fermion- loop effective action (|10|) . one then obtains 

•^frxc^l-l /"d4rrA + rT^2_.rj. , /P0-K)(x)-e[P + V(a;)] -[A + a(x)] \ , . 

where V^(a:) = d^ip{x) formally plays the role of an Abelian gauge field. In this guise, the Nambu-Goldstone field 
is invariably accompanied by at least one derivative. The resulting effective theory is thus automatically invariant 
under global U(l) phase transformations, under which tf{x) is shifted by a constant, 

(^(x) — J> 1^9(2;) 4- const. (27) 



By construction, the a field appears only in the combination A + cr. 

The leading order (LO) terms in the effective theory can be obtained by ignoring derivatives on a and so that, 
after using Eq. ([12]), ^(k±V) = ^(k)-|- V^/2m±k- V/m in Eq. (HH). We explicitly checked that the terms ±k- V/m, 
do not contribute in leading order. The constant field thus appears only in the combination X = jjL — V with 

V = Vo + ^\^ =dt^+^{V^)\ (28) 
Zra Zm 

By Eq. ([5]), the Nambu-Goldstone field transforms under a Galilei boost as 

(p{x) — ip'{x') = (p{x) — mu • X + ^mu^t, (29) 

and the two terms at the right of Eq. (|28|) combine precisely so that V is invariant. 

Given these observations, the LO terms, i.e., terms without derivatives on a and V, in the one-fermion-loop effective 
action S[a, f] = J d'^xCho + ■ ■ ■ governing these fields, must be of the form 



00 

£lo(^,¥') = E -m''^'''^'"'^' 
with TT^^J) the expansion coefficients, or vertices 

'^'-^-(-ir-^^^Ve., (31) 

which depend on the parameters to, /i, A of the theory. By dimensional analysis, 

^('.^") = _L^^3/2a5/2— (32) 

with Xo = ^J,/ A- The numerical prefactor is pulled out for later convenience. The expansion coefficients can be readily 
obtained from the expression (1131) for the effective potential. The results up to order i + j = 4 are recorded in the 
Appendix. We iterate that the only terms generated at leading order are those dictated by symmetry and are precisely 
the once included in Eq. pop . Returning to the original expression ()26p . we explicitly checked that terms spoiling any 
of the symmetries drop out. 

To the order z + j = 2, i.e., in the Gaussian approximation, the a field can be integrated out by substituting the 
corresponding Euler-Lagrange equation for this field. 
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back into the Lagrangian ([30]) with i + j <2. This is tantamount to approximating the integral over a by the saddle 
point. It gives as effective theory governing solely the Anderson-Bogoliubov mode at LO 

1 CttCI^I))^ _ (0,2) (2,0) 

£.o(^) = .(°'^V - V^. (34) 

From it, the speed of propagation c of this gapless mode can be read off as 

1 ^(0,1)7^(2.0) 

= 5 . (35 

m (7r(i,i)) _ 7r(o^2)7r(2,o) 

Substituting the explicit expressions (|A7[) for the coefficients tt*^*'^^ , we reproduce the result due to Marini, Pistolesi, 
and Strinati [ll| 

2 _ 2 /X /5(a:o/5 +/6) /ofiN, 

' 3ma;„(/| + /|)- ^^""^ 

In the weak-coupling BCS limit, obtained by letting oo, as well as in the strong-coupling unitary limit, where 

^5 = ^qIq with Xo — 0.8604 . . ., it reduces to the same form 

c' = l^- (37) 

3 TO 

Both limits are in the regime where the chemical potential is positive. At unitarity, 

1 

— = -^== 0.1968..., (38) 
w| 3 

with vp = kp /to the Fermi velocity of the free Fermi gas at the same density as the unitary gas, and ^ the dimensionless 
parameter ([22]) . 

We next turn to LO terms of higher powers in the fields. In the BCS limit Xq oo, where at leading order in 
energy and momentum, the a field decouples from V and can be ignored, the effective theory pop up to quartic order 
in V reduces to 

Com = ^,n"v" (i - f + f - |v-» - ^v-' + . . .) . m 

where we added the free Fermi gas contribution and introduced the abbreviation V = V/fi. In this limit, /i — > ep. 
The terms in Eq. p9p form precisely the first in the Taylor series expansion of the full LO expression, 

/:lo(^) = coto3/2x5/2, X = /i-y = ^-9t(^--^(V^)2, (40) 

Zm 

proposed by Son and Wingate [l^ with the coefficient 

25/2 

^° = 15^- ^41) 

The expression (j40p sums up all terms where each Nanibu-Goldstone field is accompanied by exactly one (space or 
time) derivative. For a static uniform system, where V is zero, the right side of Eq. (|40p reduces to the zero-temperature 
pressure expressed as a function of /i, P{n)- 

As an aside, a trapping or other external potential U{x) can be readily included in the formalism by replacing fi 
with II — U{x) in X. 

To obtain the effective theory governing the field V alone in the entire BCS-BEC region, the a field must be 
integrated out from the complete LO expression ([30]). Since it contains arbitrary powers of a, this is in general 
impossible. At unitary, however, the coefficients assume values that make this at least approximately possible, as we 
now demonstrate. The Lagrangian up to quartic order in a and V, with the vertices given in the Appendix, takes the 
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following form in this limit: 



Xo 



--fl + x^)V + -x^V^ - —x^V^ - - ^ 



-XoaV - ^xla^ - ^xla^ 



1 



xla^V 



1 



4 ° 



1 xl{l + 2xl) 
24 1 + a;2 



24 1 + a;2 



1 x^(3 + 4xg) 4 
96 1 + a;2 

1 : 



16 1 + j:2 



(42) 



Note that the vertex tt'^'^) of the aV^ term vanishes in this limit. Moreover, the coefficients are such that the 
Euler-Lagrange equation (|33|) . which assumes the form a — ~V/xo in the unitary limit, obtained in the Gaussian 
approximation remains unchanged after including the additional terms in Eq. (I42p . Put differently, the saddle point 
of the nonlinear theory (|42p remains locked at a — —V/xq in the unitary limit. Approximating the integral over a 
by this saddle point, we obtain as effective Lagrangian governing just ip 



(1 + 

3/2 



1 



16 



128 



(43) 



where we included the term for a static uniform unitary Fermi gas. The various contributions combine to exactly 
generate the first terms in the Taylor series expansion of the predicted form ([^0]) with 



Co 



25/2 ^ 
157r2 ^3/2 ' 



1 

^372 



3/2 



2.203. 



(44) 



so that Co = 0.0841 . . .. The dimensionless parameter ^, which was introduced in Eq. (1^^ as the ratio of the chemical 
potential and ep, gives here the ratio of the ground-state energy per particle e of the unitary gas and that of a free 
Fermi gas, or equivalently, a neutral BCS superconductor at the same density. 



c 3 kp 
^10 m' 



(45) 



as follows from taking the Legendre transform of P(/u) and using that n — dP{ii)/dn [Tjl- We emphasize that in 
determining the coefficients of the LO terms only the effective potential describing a static uniform system is 
used. This validates the symmetry argument by Son and Wingate [l^ that the LO effective theory (PH)) to all orders 
in the Nambu-Goldstone field can be obtained by simply replacing the chemical potential ^ with the Galilei-invariant 
combination X in the pressure -P(/i) of the static uniform system. It implies that the complete interaction between 
phonons at leading order in wave vector and frequency is determined by P{y)- 

The strong-coupling unitary limit is special as for no other point in the BCS-BEC crossover, the saddle point ([33|), 
obtained in the Gaussian approximation, constitutes a solution to the Euler-Lagrange equation when additional terms 
in the Lagrangian ([50)1 are included. 



IV. ID 



Although the coefficients of the LO effective Lagrangian in three space dimensions (3D) can only be determined 
approximately in the strong-coupling unitary limit, the form of the theory is precisely as predicted by Son and 
Wingate p^ . To provide further support for this prediction, we in this section consider the pairing theory in one 
space dimension (ID) for which exact results are available. 

The ID system of N spin-^ fermions with attractive delta- function interactions is described exactly by the Gaudin 
integral equations [28]. Let, as in the three-dimensional case, A(< 0) denote the coupling constant. It is related to 
the ID scattering length through 

\ = -\ma, (46) 

as follows from Eq. ^ with D = 1. Whereas the 3D coupling constant is directly proportional to the scattering length, 
its ID counterpart is inversely proportional to a, and a — > -|-cxd in the limit A — > 0~. This divergence of the scattering 
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length in the zero-couphng hmit arises because an attractive delta-function potential possesses a two-particle bound 
state however small the attraction. In other words, this limit marks the threshold of a bound state in vacuum, where 
the scattering length diverges and changes from negative for repulsive interactions (A > 0) to positive for attractive 
interactions (A < 0). 

In the weak-coupling BCS limit A 0^, the Gaudin integral equations yield for the ground-state energy per particle 
e, chemical potential fj., and sound velocity c, the free Fermi gas expressions 

e=i^, fi = eF, ^ = 2^^vl, (47) 

TO TO 

as expected. In the strong-coupling limit A — cxd, where a — > 0+ in ID, the fermions form tightly bound pairs with 
binding energy ea = 1/ma?, as in 3D. The Gaudin integral equations give in this limit [29j 

ecff - — , Meff - g — , c -2 -7 — ' (.48) 

24 TO 8 TO TO 4 TO^ 

where we removed the (diverging) binding energy from the ground-state energy and the chemical potential by intro- 
ducing Ecff = + e and /ies = + M- The Fermi momentum fcp is now defined through 

2 

n = -kp. (49) 

TT 

Note that as in 3D, the unitary limit is in the region with positive (effective) chemical potential. On comparison with 
the weak-coupling results (|T7)) . it follows that 

= J (50) 

exactly. As in 3D, ^ also equals the ratio of the (effective) ground-state energy per particle e of the unitary gas and 
that of a free Fermi gas at the same density. The exact ID counterpart of Eq. psp reads 



whereas the full LO Lagrangian density is given by 

/:lo(<p) =coto1/2x3/2, (52) 

with 

25/2 1 ^7/2 

^0 = ^^-^ (53) 

exactly in the unitary limit (and Cq = 2^/^/37r in the weak-coupling limit). 

The physical interpretation of the ID unitary limit follows from comparison with the related problem of repulsively 
interacting bosons. The ID bosonic system with a delta-function potential is described exactly by the Lieb-Liniger 
integral equations [SO^. In the infinite-coupling limit, these equations coincide with the Gaudin integral equations at 
unitary, provided one identifies the boson mass TOb with twice the fermion mass, TOb = 2to, and the boson number 
density rib with half the fermion number density, nb — n/2 29]. That is, the ID unitary limit coincides with the BEG 
limit. 

The two sets of integral equations can be mapped onto each other not only in the infinite-coupling limit, but also for 
finite coupling. With the above identifications, both sets become similar with one important distinction that the signs 
of the coupling constants differ, as already pointed out by Gaudin (28j . Specifically, the Gaudin integral equations can 
be obtained from the Lieb-Liniger integral equations by setting TOb = 2to and Tib ~ n/2, and by replacing the bosonic 
coupling constant Ab(> 0) with — 2A(< 0) showing that the interaction between pairs is attractive. Whereas 
bosons must interact repulsively to guarantee stability, pairs of spin-^ particles can have attractive interactions, for 
the Pauli exclusion principle forbids two such pairs to form a four-fermion bound state. The fact that is positive 
in the unitary limit implies that the compressibility is also positive and that the system is mechanically stable even 
in this limit of infinite attraction between pairs. 
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V. NEXT-TO-LEADING ORDER 



We next turn to the next-to-leading order (NLO) terms in the effective Lagrangian. These terms involve derivatives 
of the a field and V , or equivalently, X = ^ — V . To obtain the first NLO terms it suffice to consider quadratic terms 
in the fields cr and up to fourth order in derivatives. The possible independent terms of this form are given in 3D 
by 



,(2) 
-NLO 



1 



21/2 67^2 



A 



3/2 



1/2 



1/2 



/ fY) \ 1 /2 / jj) \ 3/2 / \ 1/2 



(54) 



where the powers of m/A follow from dimensional analysis. In writing this general expression, with arbitrary coeffi- 
cients 6i,Ci, and di, we also used that the theory is invariant under time reversal, under which t ^ —t and (p —^p. 
Our choice of coefficients is such that at unitarity, the terms with coefficients di , d2 and ^4 combine with those with 
coefficients ci,C2 and C4, respectively after the a field has been integrated out. In addition to the (V^</j)^ term, 
there exists a second term quartic in derivatives, viz. didjipdidjLp. In the quadratic approximation we are working, 
both differ by a total derivative and cannot be uniquely identified. However, the ratio of the two terms can be 
determined by the derivative expansion method we use and comes out to be two, leaving two possible combinations 
(V^tp)^ -|- 2{didj(p)^ or 2(V^(/5)^ -I- {didjipY . Through the study of two three-point correlation functions, it was shown 
in Ref. (32j that the former combination is in fact realized. This combination is also favored by symmetry, for 

{V'^ipY + 2{didjpY = {SijSki + SikSji + Sii6kj)didjpdkdi'p. (55) 

We have computed the coefficients appearing in the Lagrangian ([5l)) for the entire BCS-BEC crossover by applying 
the derivative expansion method to the formal expression (|26)) . The results are collected in the Appendix. Although 
the terms {'Vdt'p)'^ (with coefficient ci) and dfpV'^Lp (with coefficient C4) become identical after partial integration, 
they can be separately identified, for they originate from different parts in the the effective action (P^ . Specifically, 
the term {VdifpY arises as (VVg)^ and can be calculated by setting V to zero in Eq. (pS)) and ignoring further time 
derivatives, while the term dfpiV'^ip arises as dtW ■ V. 

We first consider the BCS limit, where the a field decouples in first approximation. The general expression (|5T 
reduces in this limit, obtained by letting Xo 00, to the known result (in conventional notation) [TgI, [2ll. [ssj 



(2) 
NLO 



if) 



6 A2 



(dt^y + iv^v^py^^-v^dtpv^p 



(56) 



with i^{0) = mkp/2TT'^ the density of states at the Fermi surface. Note that in this limit, the coefficient ci can be 
ignored in comparison to the coefficient C4. These NLO terms together with the quadratic LO terms in Eq. ([M)) give 
rise to the phonon spectrum [l^, [sl] 



1 



2i 2 



1 



^ '^F 1 2 

45 A2 



2 ^, 



3 V 45 A2 J 3 m 
The quadratic terms (j56l) can be put into Galilei-invariant form as 



^NLo('y5) = Y^coxl 



,3/2 



:{DtXf 



4 

15 ml/2 



^^k2 

45 TO/i 



4 ml/2 
3X1/2 



DtXV^p 



(57) 



(58) 



with Co given in Eq. (|41|) and Dt the material derivative. 



Dt = dt + -V(^.V , 
m 



(59) 



which by the transformations ([4]) and ((29|) is invariant under Galilei boosts. We note that the coefficients of the NLO 
terms can be uniquely identified from the quadratic approximation (I56p to this Lagrangian. Indeed, as prescribed by 
the derivative expansion method [24[, each factor of jjl is to be replaced with X in Eq. (|58l) . And each occurrence of 
dfip is to be replaced with the Galilei-invariant form —DtX. 



10 



For completeness, we mention that the Lagrangian governing the a field in the BCS limit reads to this order 



1-2 1-3 
(T rCT 

2 3! 



4! 



(60) 



where a = cr/A, and vanishes in the limit Xo — >■ oo. 

We continue by studying the strong-coupling unitary limit obtained by setting Ic^ = XqIq with Xo = 0.8604 . . .. The 
general expression (I54p assumes in this limit the explicit form: 



'^NLo(^'</') = 



1 



21/23^2 



1 ( m\^^^ \2 I I m 

4 7 -oid.ar-T4^- 



1/2 



1 



1/2 



1/2 



2"° 



(7 + 4x2)(Va)2 
3/2 



- ) (d^^r 



+ |(l + -o) (^)'^'(VW- 1(4x^ + 1) 
5 Vto/ 6 



1/2 



3/2 



1/2 



XodtcrV'^Lp 



(61) 



Note that the vertex di of the Vcr • Vdtip term vanishes in this limit. Approximating the integral over a by the LO 
saddle point p3p . which is consistent to the order we are working, we obtain an effective theory of precisely the form 
(j56p with an additional {Vdt^p)^ term included. As in the BCS limit, these quadratic terms fourth order in derivatives 
can be put into Galilei-invariant form as 



7 1 ml/2 
12^X1/2 



{VXf 



,3/2 



X3/2 



4 X^/^ 
15 m-'^/^ 



[(vV)' + md,'^? 



4mi^ 
3XV2 



DtXV^Lp 



(62) 

with Co now given by Eq. (j44p . Apart from overall normalization and the first term, which is suppressed in the BCS 
limit, these NLO terms are exactly as found in the BCS limit. As in that limit, the form and coefficients of the NLO 
terms uniquely follow from the quadratic approximation to this Lagrangian. In addition to the replacements already 
used in the BCS limit, (Vdi^pY replaced with the Galilei-invariant form (VX)^ in Eq. The second and last 

terms in that expression, both involving DtX, were not considered by Son and Wingate [13] as their approach is 
limited to only gradients of X. Also the fourth term was omitted. With {i — 1,2,3,4,5) denoting the coefficients 
of the NLO terms, so that 



7 1^4 



8 4 



:Coa;; 1 - — — , i, — , — , 77 ) , (63) 
the spectrum of the gapless Anderson-Bogoliubov mode that follows when including the NLO terms reads 



2/1 \ 2 /i 2 

"('')=3m^ 

^2^j^2 

3 m 



1 



45co 

35 - 32x2 



[6ri + 4r2 + 9{r3 + n) 
1 , 



360 m/i 



(64a) 
(64b) 



It reduces to the BCS expression (|57)) in the limit Xo 00. In contrast to the BCS limit, the coefficient of the 
correction term is positive for the value Xq = 0.8604. . . obtained in the saddle-point approximation. If it remains 
positive beyond this approximation, a low-energy phonon in a unitary Fermi gas can decay into two phonons. 

The Son- Wingate result for the spectrum corresponds to setting r2 = = and also r4 = in Eq. (|64ap . which 
leads to an incorrect expression for the spectrum. The correct result (j64bp was derived from a NLO Lagrangian of a 
form proposed by Son and Wingate, i.e., one without time derivatives in Ref. [32*1 . This was achieved by eliminating 
the time derivatives in the NLO Lagrangian through the use of the leading-order field equations. This reduction leads 
to the following changes in the coefficients (|63)) : 



1 2 
7-3 = Y^coa;„ -> 



1 2 

'18^°"°' 



r4 = 2r3 



= -3r' 



(65) 



which, when substituted in Eq. (|64ap with r2 = ^5 = 0, yields the correct spectrum (j64bp . The flip side of this 
reduction is that the static response functions come out incorrectly, for the coefficients of the static terms in the NLO 
Lagrangian now also include dynamic effects. The relation r'^ = —3r'^ in the reduced NLO Lagrangian was argued in 
Ref. 32] to be a consequence of conformal invariance at unitarity. 



11 



VI. DISCUSSION 



Its (relatively) simple form and the very fact that the effective action of a unitary Fermi gas up to next-to-leading 
order can be derived from the microscopic theory analytically underscores the special status of the unitary limit in 
the BCS-BEC crossover. Although the coefficients could only be computed approximately, using an (extended) saddle 
point, it is remarkable that the effective field program, which involves integrating out the fermionic degrees of freedom 
as well as the a field, can be carried out consistently up to the orders considered, featuring higher-order terms such 
as (Vi/?)® and {d^(p)'^. Surprisingly, the effective actions obtained in the weak-coupling BCS and the strong-coupling 
unitary limits are proportional, save for a next-to-leading order term which is suppressed in the BCS limit. 
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Appendix A: Integrals 



All integrals encountered in this study can be expressed as linear combinations of two basic integrals introduced by 
Marini, Pistolesi, and Strinati [ill 



h{xo) 



dx 



(Al) 



by integrating by parts and simple algebraic manipulations. Here, x,^x, and denote the dimensionless variables 



2 _ P/2m ^ _ ( 



A 



X = — E ^ — - JP + 1 

o, — ^ , j^x — ^ ~ V I 



(A2) 



These i nteg rals can be expressed in terms of the complete elliptic integrals of the first and second kind, as was done 
in Ref. [llf or in terms of Legendre functions Pa [1^ as (70 = Xoj \J\ -|- x^) 



4 (l-Ha;2)3/4 

TT 1 



4 (1-f x2)3/4 



(1 - 3x^)Fi/2 (-70) - 3.T<,(1 + xlf'^P^,2 {~lo) 
4XoPi/2 (-7o) + 3(1 + xlY'^Py2 (-7o) 



(A3) 



so that, for example. 



h{Xo) - XohiXo) = ^(1 + xiy/^Pi/2 (-7o) 



(A4) 



With the definition 



d-^^ (A5) 
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so that /s = /2,o,3 and Iq = /2.1.3, one readily verifies the relations 

^0,04 

^0,0,3 



2/6 



/n.n.5 — 



4, 



1 + 2:2 
-^4,0,5 ^ 2/24,5 + /o,0,3 

/s - So/e 



1,3 



2,5 



1 + 2:2 
-^0,0,3 ~ -^0,0,5 
-^0,3,5 = -^2,2,5 — 2:0/0,2,5 

-^2,0,5 



1 2:0/5 + -^6 , 1 ^ 



^2,0,7 — ^2,4,7 + 2/2,0,5 — -^5 
1 



^2,1,5 
-^2,1,7 
^2,2,5 
^2,4,7 
^4,0,5 



:/( 



^-0,0,3 

6 

^/0,0,5 
0^5 + 77^0,1,3 



1 

10 



6 

(-^^0,3,5 + 6/2,2,5) 
1 



(1 + a:o)/2,i,5 + -2:0/5- 



These integrals are all a function of Xo alone. 

With the help of these integrals, the coefficients tt^*'-'^ of the LO effective theory (|30|) . introduced in Eq. 
be readily evaluated in closed form, with the results up to i + j = 4 



,(0,1 



r(0,2 



f(0,3 



f(0,4 



f(3,0 



f(4,0: 



f(1.2 



.(1,3 



f(2,l 



r(2,2 



,(3,1 



2/5 
20/5 + h 

1 (-3 + xl)h + 4x0/6 

2 (l + ^?)2 

-2/5 

(3 + 2a:g)/5 - Xq/q 
1 + 2:2 

1 (3 + 32:g + 42:^)4 + 4a:g/6 

2 (1 + 2:2)2 

-2/6 

h - Xoh 

1 + xl 

1 4x0/5 + (1 - 32:g)/6 

2 (l + a:2)2 

X0/5 - (1 + 2xg)/6 

l(l-3xg)/5+2xo(-l + xg)/6 

2 (1 + ^:2)2 

1 2xo(-l + xg)/5 + (l + 5xg)/6 
2 (l + a:2)2 



(A6a) 
(A6b) 

(A6c) 

(A6d) 

(A6e) 
(A6f) 

(A6g) 

(A6h) 

(A6i) 

(A6j) 

(A6k) 

(A61) 

(A6m) 

32|) . can 

(A7a) 
(A7b) 
(A7c) 

(A7d) 
(A7e) 
(A7f) 

(A7g) 
(A7h) 
(A7i) 

(A7j) 

(A7k) 

(A71) 

(A7m) 
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Note that 7f(0'2) ^ _^(2,o)_ 

The coefficients of the NLO terms appearing in the quadratic Lagrangian ([5^ . which are somewhat laborious to 
compute, can again be expressed as linear combinations of and /g, with the results 

, 1 a:„(l+4a;g)/5 + (7 + 10xg)/6 

- "12 TT^l (A8b) 

1 X0I5 +h , , „ ^ 

^^^l(3 + 4x|)4+x^ (A8d) 



4 1 

^1 



c:i = \[{i + ^xl)h+Xoh] (A8e) 



c^^-^(Axoh+h) (A8f) 
di = I^. (A8i) 
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